Squiral Diffraction 



Uwe Grimm and Michael Baake 



Abstract The Thue-Morse system is a paradigm of singular continuous diffraction 
in one dimension. Here, we consider a planar system, constructed by a bijective 
block substitution rule, which is locally equivalent to the squiral inflation rule. For 
balanced weights, its diffraction is purely singular continuous. The diffraction mea- 
sure is a two-dimensional Riesz product that can be calculated explicitly. 



1 Introduction 

The diffraction of crystallographic structures and of aperiodic structures based on 
cut and project sets (or model sets) is well understood; see J4] [5] and references 
therein. These systems (in the case of model sets under suitable assumptions on the 
window) are pure point diffractive, and the diffraction can be calculated explicitly. 

The picture changes for structures with continuous diffraction. Not much is 
known in general, in particular for the case of singular continuous diffraction, even 
though both absolutely and singular continuous diffraction show up in real systems 
US E). The paradigm of singular continuous diffraction is the Thue-Morse chain, 
which in its balanced form (constructed via the primitive inflation rule 1 i — >- 11, 
1 i — y 1 1 with with weights 1 and I = — 1, say) shows purely singular continuous 
diffraction. This was shown by Kakutani iflOl . see also fTJ, and the result can be 
extended to an entire family of generalised Thue-Morse sequences J3|. 

Here, we describe a two-dimensional system which, in its balanced form, has 
purely singular continuous diffraction. For more detail and mathematical proofs, we 
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refer to J6). Again, it is possible to obtain an explicit formula for the diffraction 
measure in terms of a Riesz product, with convergence in the vague topology. 



2 The squiral block inflation 

The squiral tiling (a name that comprises 'square' and 'spiral') was introduced in J8] 
Fig. 10.1.4] as an example of an inflation tiling with prototiles comprising infinitely 
many edges. The inflation rule is shown in Fig. [T] it is compatible with reflection 
symmetry, so that the reflected prototile is inflated accordingly. 



Fig. 1 The primitive inflation 
rule for the squiral tiling of 
the Euclidean plane. 



Jl 



A patch of the tiling is shown in Fig. [2] Clearly, the tiling consists of a two- 
colouring of the square lattice, with each square comprising four squiral tiles of the 
same chirality. The two-colouring can be obtained by the simple block inflation rule 
sown in Fig. [3] which is bijective in the sense of ifTTl . Again, the rule is compatible 
with colour exchange. The corresponding hull has D\ symmetry, and also contains 
an element with exact individual D4 symmetry; see |6J for details and an illustration. 

Due to the dihedral symmetry of the inflation tiling, it suffices to consider a tiling 
of the positive quadrant. Using the lower left point of the square as the reference 
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Fig. 2 Patch of the squiral tiling obtained by two inflation steps from the central seed. 
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Fig. 3 Equivalent block infla- 
tion rule for the squiral tiling 
of Figure [2] 



□ 



□ 



point, the induced block inflation p produces a two-cycle of configurations v and 
pv. They satisfy, for all m,n > and < r,s < 2, the fixed point equations 

_ ifr = s = 0mod2, 

{Pv)3m+r3n+s - ) v otherwise. (l) 



3 Autocorrelation and diffraction measure 

For a fixed point tiling under p 2 , we mark each (coloured) square by a point at its 
lower left corner z <G Z 2 . For the balanced version, each point carries a weight w T = 1 
(for white) or w- = 1 = — 1 (for grey). Consider the corresponding Dirac comb 



' S Z 2 = H W z S z- ( 2 ) 



CO 



Following the approach pioneered by Hof |5), the natural autocorrelation measure 7 
of o is defined as 

r -..»^ fe W|^Wa) , (3 ) 

where Cn stands for the closed centred square of side length 2iV. Here, denotes 
the measure defined by fi.(g) = forg 6 C C (R 2 ), with := g(— x) (and where 
the bar denotes complex conjugation). The autocorrelation measure 7 is of the form 
7= 775^2 with autocorrelation coefficients 

1 N 

7} (m, n) = ton y N+1 y ^ L W U W k-m,t-n ■ ( 4 ) 

All limits exists due to the unique ergodicity of the underlying dynamical system 
|6|, under the action of the group Z 2 . 

Clearly, one has 77(0,0) = 1, while Eq. ([TJ implies the nine recursion relations 

77(3m,3n) = T](m,«), 
77(3m,3n+l) = -lr\{m,n) + \r\(m,n+\) , 
Tj(3m,3n+2) = \r\(m,n) - %T}(m,n+l) , 
77(3m+l,3«) = -^r\(m,n) + \r]{m+l,n), 
Tj(3m+1,3«+1) = -§(»](m+l,n) + T/(m,n+l)) + ^(m+l,n+l), (5) 
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rj(3m+l,3n+2) = -|(Tj(m,n) + Tj(m+l,n+l)) + ^T}(m+l,n), 

rj(3m+2,3n) = \r\{m,n) - |Tj(m+l,n) , 

rj(3m+2,3n+l) = -|(Tj(m,n) + Tj(m+l,n+l)) + iTj(m,n+l), 

t7(3m+2,3«+2) = \r\{m,n) - %(r\(m+\,n) + r\{m,n+l)), 

which hold for all m,n £ Z and determine all coefficients uniquely (6|. The autocor- 
relation coefficients show a number of remarkable properties, which are interesting 
in their own right, and useful for explicit calculations. 

Since the support of (0 is the lattice Z 2 , the diffraction measure y is Z 2 -periodic 
Q, and can thus be written as 

r = n * <s Z 2 , 

where fx is a positive measure on the fundamental domain T 2 = [0, l) 2 of 1? . One 
can now analyse y via the measure /i, which, via the Herglotz-Bochner theorem, is 
related to the autocorrelation coefficients by Fourier transform 

TJ(*)= I / K,kz dn{z), 

Jt 2 

where k = (m,n) & 1? and kz denotes the scalar product. We now sketch how to 
determine the spectral type of jj., and how to calculate it. 

Defining £(Af) := ^y^ 77 (m,n) 2 , the recursions © lead to the estimate 

E(3N)<^-E(N), 

so that E(N)/N 2 — > as —> °°. An application of Wiener's criterion in its multi- 
dimensional version (6][7) implies that jj,, and hence also the diffraction measure y, 
is continuous, which means that it comprises no Bragg peaks at all. 

Since I] (0,1) = 17 (1,0) = —1/3, which follows from Eq. (O by a short calcu- 
lation, the first recurrence relation implies that 77 (0, 3 J ) = Tj(3', 0) = — 1/3 for all 
integer j > 0. Consequently, the coefficients cannot vanish at infinity. Due to the 
linearity of the recursion relations, the Riemann-Lebesgue lemma implies |6| that 
jj. cannot have an absolutely continuous component (relative to Lebesgue measure). 
The measure jj., and hence y as well, must thus be purely singular continuous. 



4 Riesz product representation 

Although the determination of the spectral type of y is based on an abstract argu- 
ment, the recursion relations (O hold the key to an explicit, iterative calculation of 
/X (and hence y). One defines the distribution function F(x,y) := ji([0,x] x [0,y]) 
for rectangles with < x,y < 1, which is then extended to the positive quadrant as 

F(x,y) = y([0,x]x[0,y]). 
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This can finally be extended to R 2 via F(—x,y) = F(x, —y) = —F(x,y) and hence 
F(— x,— y) =F(x,y). In particular, one hasF(0,0) =0as well a&F(Q,y) =F(x,Q) = 
0, and F is continuous on R 2 . The latter property is non-trivial, and follows from 
the continuity of certain marginals; see J6) and references therein for details. 
One can show that, as a result of Eq. (0, F satisfies the functional relation 

written in Lebesgue-Stieltjes notation, with the trigonometric kernel function 

1 2 

T^(jc,y) = — (l +2cos(27Dc) +2cos(27ry) — 4cos(27Dc) cos(27Ty)) . 

The functional relation © induces an iterative approximation of F as follows. Start- 
ing from F^(x,y) = xy (which corresponds to Lebesgue measure, dF^ = X) and 
continuing with the iteration 

F {N+l) ( X ,y) = i f X [ 3y H^)dFW(x,y), (7) 
9 Jo Jo 3 3 

one obtains a uniformly (but not absolutely) converging sequence of distribu- 
tion functions, each of which represents an absolutely continuous measure. With 

dFM (x,y) = f (N) (x,y) dxdy, where /W (x,y) = -^F (N) (x,y), one finds the Riesz 
product 

f iN) (x,y) = Y\&(3 e x,3 e y). (8) 

These functions of increasing 'spikiness' represent a sequence of (absolutely contin- 
uous) measures that converge to the singular continuous squiral diffraction measure 
in the vague topology. The case N = 3 is illustrated in Figure|4] Local scaling prop- 
erties can be derived from Eq. ((H). 
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The example of the squiral tiling demonstrates that the constructive approach of 
Refs. iflOl [U O can be extended to more than one dimension. The result is as ex- 
pected, and analogous arguments apply to a large class of binary block substitutions 
that are bijective in the sense of ifTTI . This leads to a better understanding of binary 
systems with purely singular continuous diffraction. 

It is desirable to extend this type of analysis to substitution systems with larger al- 
phabets. Although the basic theory is developed in |Q~2), there is a lack of concretely 
worked-out examples. Moreover, there are various open questions in this direction, 
including the (non-)existence of bijective constant-length substitutions with abso- 
lutely continuous spectrum (the celebrated example from lfl"2l Ex. 9.3] was recently 
recognised to be inconclusive by Alan Bartlett and Boris Solomyak). 
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